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We present a theoretical study using density functional calculations of the structural, electronic 
and magnetic properties of 3d transition metal, noble metal and Zn atoms interacting with carbon 
monovacancies in graphene. We pay special attention to the electronic and magnetic properties of 
these substitutional impurities and found that they can be fully understood using a simple model 
based on the hybridization between the states of the metal atom, particularly the d shell, and the 
defect levels associated with an unreconstructed Ds^ carbon vacancy. We identify three different 
regimes associated with the occupation of different carbon-metal hybridized electronic levels: (i) 
bonding states are completely filled for Sc and Ti, and these impurities are non-magnetic; (ii) the 
non-bonding d shell is partially occupied for V, Cr and Mn and, correspondingly, these impurties 
present large and localized spin moments; (Hi) antibonding states with increasing carbon character 
are progressively filled for Co, Ni, the noble metals and Zn. The spin moments of these impurities 
oscillate between and 1 and are increasingly delocalized. The substitutional Zn suffers a Jahn- 
Teller-like distortion from the C^v symmetry and, as a consequence, has a zero spin moment. Fe 
occupies a distinct position at the border between regimes (ii) and (Hi) and shows a more complex 
behavior: while is non-magnetic at the level of GGA calculations, its spin moment can be switched 
on using GGA-f-U calculations with moderate values of the U parameter. 

PACS numbers: 73.22.-f, 73.20.Hb, 75.20.Hr, 61.48.De 



I. INTRODUCTION 



The electronic properties of two-dimensional graphene 
are currently the subject of intense experimental and 
theoretical researchi^^ Graphene exhibits many intrigu- 
ing phenomena stemming from the characteristic conical 
dispersion and chiral behavior of its valence and con- 
duction bands nearby the Fermi energy^^i^i^ In general, 
graphenic nanostructures like graphene nanoribbons, car- 
bon nanotubes and their interconnections, are opening 
new routes for research in the field of nanoelectronics^, 
due to the long spin relaxation and decoherence times of 
these materials^. In addition, the possibility to control 
the magnetism of edge states in nanoribbons by apply- 
ing external electric fields^i^ make them very attractive 
for spintronic applications^^ However, for the design of 
realistic devices the effect of defects and impurities have 
to be taken into account. For this reason, quite a lot 
of work has been devoted to the study of defects and 
different types of impurities in these materials. It has 
now become clear that defects and dopants severely af- 
fect some the properties of graphenic systems and can be 
used to tune their response. For example, the strong in- 
terplay between the presence of defect or dopants and the 
magnetic properties of carbon nanostructures has been 
stressed by several authors piiii^ii^iiiii^ii^iiiii^ii^ 

Here we focus on substitutional impurities in graphene, 
in which a single metal atom substitutes one or several 
carbon atoms in the layer. Direct experimental evidence 
of the existence of these kind of defects has been recently 
provided by Gan et al^. Using high-resolution transmis- 
sion electron microscopy (HRTEM), these authors were 



able to visualize individual Au and Ft atoms incorpo- 
rated into a very thin graphitic layer probably consist- 
ing of one or two graphene layers. From the real-time 
evolution and temperature dependence of the dynamics 
they obtained information about the diffusion of these 
atoms. Large diffusion barriers (~2.5 eV) were observed 
for in-plane migration, which indicates the large stability 
of these defects and the presence of strong carbon-metal 
bonds. These observations seem to indicated that the 
atoms occupy substitutional positions. There also ex- 
ists evidence for the existence of substitutional Ni impu- 
rities in single-walled carbon nanotubes (SWCNT)2ii^ 
and graphitic particlesj^ Ushiro et al^ showed that Ni 
substitutional defects were present in SWCNT samples 
even after careful purification and that, according to their 
analysis of the X-ray absorption data, the most likely 
configuration was one in which the Ni atom replaces a 
carbon atom, as seen in Fig. [T] 

The presence of substitutional defects can have impor- 
tant implications for the interpretation of some experi- 
mental evidence. For example, substitutionals of mag- 
netic transition metals are expected to strongly influence 
the magnetic properties of graphenic nanostructures. In- 
terestingly transition metals, like Fe, Ni or Co, are among 
the most common catalysts used for the production of 
SWCNT^"'. Furthermore, Rodriguez-Manzo and Ban- 
hart^^ have recently demonstrated the possibility to cre- 
ate individual vacancies at desired locations in carbon 
nanotubes using electron beams. This ability, in combi- 
nation with the observed stability of substitutional im- 
purities, can open a route to fabricate new devices in- 
corporating substitutional impurities in certain locations 
or arranged in particular ways. Such devices would al- 
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low for experimental verification of some of the unsual 
magnetic interactions mediated by the graphenic carbon 
network that have been predicted recently.— 

In this work we study theoretically the structural, 
electronic and magnetic properties of 3d transition met- 
als, noble metals and Zn as substitutional impurities in 
graphene using density functional (DFT) calculations. 
We only consider the configuration proposed by Ushiro 
et al.,^^ in which a single metal atom binds to a car- 
bon monovacancy. Throughout the paper we will refer 
to this structure as substitutional configuration. The 
results of our DFT calculations are in good agreement 
with the results of other recent studies on similar sys- 
tems j^'^^'^S'^i^'^ In particular, Krasheninnikov et al^ 
have presented the most complete first-principles study 
to date of metal atoms interacting with single and double 
vacancies in graphene. In most cases our results are in 
good agreement with their predictions and similar gen- 
eral trends are found. However, the electronic structure 
of these defects has not been analyzed in detail to date 
and no simple model to understand the observed behav- 
iors has been presented so far. Here we present the evo- 
lution of the electronic structure of the substitutional de- 
fects as we move along the transition series and correlate 
the observed changes with a simple picture of the metal- 
vacancy interaction. 

In our study we pay special attention to the electronic 
and magnetic properties of the metal impurities. One 
of our key results is that the electronic and magnetic 
properties of these substitutional metals can be fully un- 
derstood using a simple model based on the hybridiza- 
tion between the states of the metal atom, particularly 
the d shell, and the defect levels associated with the un- 
reconstructed carbon vacancy. The predictions of this 
model are in good agreement with the calculated DFT 
band structures. With this model we can easily under- 
stand the non-trivial behavior found for the binding en- 
ergy and for the size and localization of the spin mo- 
ment as we increase the number of valence electrons in 
the impurity. In brief, we have identified three different 
regimes that can be correlated with the electron filling 
of different carbon- metal hybridized levels: (i) bonding 
states are filled for Sc and Ti, and these impurities are 
non-magnetic; (ii) the non-bonding d shell is partially 
occupied for V, Cr and Mn and, correspondingly, these 
impurities present large and localized spin moments; [iii) 
antibonding states with increasing carbon character are 
progressively filled for Co, Ni, the noble metals and Zn, 
giving rise to spin moments that oscillate between and 
1 fiB and are increasingly dclocalized. 

We have also found that Zn becomes non-magnetic due 
to a Jahn- Teller distortion. However, it is possible to 
stabilize a symmetric configuration with a spin moment 
of 2fiB with a very small energy penalty of ~ 150 meV. 
Finally, our calculations confirm that, the unexpected re- 
sult that Au substitutionals^^''^*' present a spin moment of 
also holds for Ag and Cu, and thus stems only from 
the number of valence electrons (see the Slater-Pauling- 



like plot in Fig. [3]). 

We have studied with special detail the complex case 
of Fe. This impurity occupies a distinct position at the 
boundary between two different regimes and its mag- 
netic behavior stems from the competition between the 
carbon-metal hybridization and the electron-electron in- 
teraction within the 3d shell. As a result, the spin mo- 
ment of Fe is specially difficult to describe: although non- 
magnetic using standard functionals within the general- 
ized gradient approximation (GGA)^'^, the magnetism of 
Fe appears using the so-called LDA-I-U methodology with 
reasonably low values of the U parameter. 

The paper is organized as follows. After a brief de- 
scription of the computational approach in Sec. |TT1 we 
present a summary of the structure, energetics and mag- 
netic properties of all the studied elements in Sec. lIIII In 
this section we also indicate the general ideas behind our 
model of the metal-carbon hybrization in these systems. 
In Sec.|IV]the electronic structure of the unreconstructed 
D3/1 carbon vacancy in graphene is presented. This is one 
of the key ingredients to understand the binding and elec- 
tronic structure of substitutional impurities in graphene. 
The electronic structure of the different groups of impu- 
rities is described in the subsequent sections, particularly 
in Sec. [V] The Zn substitutional impurity with its Janh- 
Teller distorsions is described in Sec. IVII A special sec- 
tion (Sec. IVIip is devoted to describe the special role of 
Fe at the border between two different regimes. Finally, 
we close with some general conclusions. 



II. THEORETICAL METHODS 

We have used two different approaches to perform 
our DFT calculations: the Siesta method^-i^^i^ us- 
ing a basis set of localized numerical-atomic-orbitals^Si^ 
(NAOs) and the Vasp code^i^ using a basis set of plane 
waves. We have used the Perdew-Burke-Ernzerho^ 
GGA (PBE-GGA) functional in all our calculations. 
Most of our results are obtained using a 4x4 supercell. 
This supercell is sufficiently large to obtain reliable re- 
sults. This has been proven by performing calculations 
using larger supercells up to 8x8 for several elements. 
The utilized codes perform three-dimensional periodic 
calculations. In order to avoid spurious interactions be- 
tween periodic images of the defective graphene layer, 
the size of the supercell perpendicular to the plane was 
always larger than 15 A. The convergence with respect 
to the number of fc-points was specially critical to ob- 
tain accurate results for the spin moment in the studied 
systems. For all the impurities and different supercell 
sizes we have used a large number of fc-points, consistent 
with a 136x136 Monkhorst-Pack^° sampling of the unit 
cell of graphene, in combination with a Fermi-smearing 
of 21 meV. All the atomic coordinates were always opti- 
mized until forces along all directions were smaller than 
0.05 eV/A. 

For the Siesta calculations we have used TrouUier- 
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Martins norm-conserving pseudopotentials - generated 
using the pseudization radii shown in the Appendix. The 
pseudopotentials for the metal atoms include nonlinear 
core corrections''^ for exchange and correlation. The 
pseudocore radii {vcore) have been optimized for each 
element and are also presented in the Appedix. Using 
nonlinear core corrections is known to be critical to de- 
scribe properly the spin moment and magnetic properties 
of transition metals. We have tested that these pseu- 
dopotentials yield the correct spin moments and band 
structures in bulk phases. 

The spacing of the real-space grid used to calculate 
the Hartree and exchange-correlation contribution to the 
total energy and Hamiltonian with Siesta was equiva- 
lent to a 180 Ry plane-wave cutoff. A double-C polarized 
(DZP)^i2i basis set has been used for the calculation 
of the magnetic and electronic properties. However, we 
have checked that using a double-C (DZ) basis set for car- 
bon yields almost identical relaxed structures as the DZP 
basis and, therefore, we have used the smaller DZ basis 
for most of the structural relaxations. The shape of the 
basis orbitals was automatically determined by Siesta 
using the algorithms described in Ref. [35j. The cutoff 
radii of the different orbitals was obtained using an en- 
ergy shift of 50 meV. Although this basis set proved to 
be sufficiently accurate to describe the geometries, spin 
moments and band structures, for some metal atoms the 
binding energies where slightly overestimated. For those 
atoms the radii of the basis orbitals were enlarged (using 
smaller values of the energy shift parameter) until bind- 
ing energies were converged within a few tens of meV. 
The resulting radii are shown in the Appendix. 

With the Vasp code we have used a well converged 
plane-wave cutoff energy of 400 eV combined with the 
projected-augmented-wave (PAW) method. Using the 
provided PAW potentials allows to check the possible 
limitations of pseudopotential calculations. However, as 
shown below, the agreement between both sets of cal- 
culations is excellent. We have also performed GGA-fU 
calculations using the formulation of Dudarev et al*^ . In 
this formulation a single U parameter is used, which we 
have taken as an empirical parameter and varied in the 
range 1-4.5 eV in our calculations for all the studied im- 
purities. Only for Fe, GGA-I-U results showed significant 
(qualitative) differences respect to those of PBE-GGA 
calculations. 



III. MAIN PROPERTIES OF 
SUBSTITUTIONAL TRANSITION METALS IN 
GRAPHENE 

We provide in this section a brief summary of our re- 
sults for the structure, binding, and spin moments of sub- 
titutional 3d transition metals, noble metals and Zn in 
graphene. 




FIG. 1: (Color online) Typical geometry of transition and no- 
ble substitutional metal atoms in graphene. The metal atom 
moves upwards from the layer and occupies, in most cases, 
an almost perfectly symmetric three-fold position with Cs^ 
symmetry. 



A. Geometry and structural parameters 

The typical geometry of the systems studied in this 
paper is presented in Fig. [TJ The metal atom appears 
always displaced from the carbon layer. The height over 
the plane defined by its three nearest carbon neighbors is 
in the range 1.7-0.9 A. These three carbon atoms are also 
displaced over the average position of the graphene layer 
by 0.3-0.5 A. The total height (h^) of the metal atom 
over the graphene plane is the sum of these two contri- 
butions and ranges between 1.2-1.8 A, as shown in panel 
(c) of Fig. [21 In most cases the metal atom occupies an 
almost perfect symmetric configuration with C^y sym- 
metry. Exceptions are the noble metals, that are slightly 
displaced from the central position, and Zn that suffers 
a Jahn- Teller distortion in its most stable configuration. 
However, we have found that it is also possible to sta- 
bilize a symmetric configuration for Zn with a binding 
energy only ~150 meV smaller. This configuration was 
overlooked in a recently published study on these sys- 
temsSi and we will refer to it as Znc3„ throughout the 
paper. 

Figures [2] presents a summary of the structural param- 
eters of substitutional transition metals, noble metals 
and Zn in graphene. Solid circles correspond to calcula- 
tions using the Siesta code with pseudopotentials and 
a basis set of atomic orbitals, while open squares stand 
for Vasp calculations using plane- waves and PAW poten- 
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FIG. 2: Structural parameters and binding energies of sub- 
stitutional transition and noble metals in graphene. Bond 
lengths and angles have been averaged for the noble metals. 
The data presented for Zn correspond to the high-spin solu- 
tion with Csv symmetry, and are very close to the averaged 
results for the most stable distorted solution. 



tials. The agreement between both sets of calculations is 
quite remarkable. Data in these figures correspond to cal- 
culations using a 4x4 supercell of graphene. For several 
metals we have also performed calculations using a larger 
8x8 supercell and find almost identical results. This is 
particularly true for the total spin moments, which are 
less dependent on the size of the supercell, but require 
a sufficiently dense k-point sampling to converge. In the 
following we will mainly discuss the results obtained with 
the smaller cell since the plots of the band structures 
are easier to interpret in that case. Finally, as already 
mentioned, noble metals and Zn present a distorted con- 
figuration. A detailed description of the structural pa- 
rameters in these cases will be given below, here we only 
present the averaged structural data for noble metals and 
those corresponding to the Znc3„ case. 

The data in Fig. [5] are basically consistent with those 
reported by Krasheninnikov et al. in Ref . [29|. The behav- 
ior of the metal-carbon bond length and the height (h^) 
of the impurity over the layer reflects approximately the 
size of the metal atom. For transition metals these dis- 
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FIG. 3: (Color online) Spin moment of substitutional transi- 
tion and noble metals in graphene as a function of the number 
of valence electrons (Slater-Pauling- type plot). Black symbols 
correspond to the most stable configurations using GGA. Re- 
sults are almost identical using Siesta and Vasp codes. Three 
main regimes are found as explained in detail in the text: (i) 
filling of the metal-carbon bonding states gives rise to the non- 
magnetic behavior of Ti and Sc; (li) non-bonding d states are 
filled for V, Cr and Mn giving rise to high spin moments; [Hi) 
for Fe all non-bonding levels are occupied and metal-carbon 
antibonding states start to be filled giving rise to the observed 
oscillatory behavior for Co, Ni, Cu and Zn. Open and gray 
(red online) symbols correspond, respectively, to calculations 
of Fe using GGA-I-U and artificially increasing the height of 
the metal atom over the graphene layer (see the text). Symbol 
marked as Zucg^ corresponds to a Zn impurity in a high-spin 
symmetric Csu configuration. 



tances decrease as we increase the atomic number, with 
a small discontinuity when going from Mn to Fe. The 
carbon-metal bond length reaches its minimum for Fe 
(dc-Fe=l-76 A), keeping a very similar value for Co and 
Ni. For Cu and Zn the distances increase reflecting the 
fully occupied 3d shell and the large size of the 4s orbitals. 
Among the noble metals we find that, as expected, the 
bond length largely increases for Ag with respect to Cu, 
but slightly descreases when going from Ag to Au. The 
latter behavior can be understood from the compression 
of the 6s shell due to scalar relativistic effects. 



B. Binding energies 

The binding energies of the studied substitutional 
metal atoms in graphene can be found in panel (d) of 
Fig. [21 In general, the behavior of the binding energies 
can be correlated with that of the carbon-metal bond 
length, although the former is somewhat more compli- 
cated. Binding energies for transition metals are in the 
range of 8-6 eV. Ti presents the maximum binding en- 
ergy, which can be easily understood since for this ele- 
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FIG. 4: (Color online) (a) Calculated band structure for an 
unrelaxed carbon vacancy (Dsh symmetry) in a 4x4 super- 
cell of graphene. Symbols indicated those bands with larger 
weight on the carbon atoms around the vacancy (solid sym- 
bols for bands with sp character and open symbols for a band 
with pz character). The electronic structure near Ef is dom- 
inated by a fully symmetric pz level (A pz) and two defect 
levels with E symmetry and sp character (E sp). Notice that, 
due to the strong hybridization with the rest of the graphene 
layer, it is not possible to identify well defined defect levels 
with E symmetry and Pz character, (b) Scheme of the elec- 
tronic structure of the Da^ C vacancy indicating the character 
and symmetry of the different levels and their occupations. 
Signal associated with the E pz level extends over the whole 
valence and conduction band. 



merit all the metal-carbon bonding states (see Sec. IV A[) 
become fully occupied. One could expect a continuous 
decrease of the binding energy as we move away from 
Ti along the transition metal series and first the non- 
bonding Sd and later the metal-carbon antibonding lev- 
els become populated. However, the behavior is non- 
monotonic and the smaller binding energies among the 
3(i transition metals are found for Cr and Mn, while a 
local maximum is observed for Co. This complex behav- 
ior is related to the simultaneous energy down-shift and 
compression of the 3d shell of the metal as we increase 
the atomic number. This will become more transparent 
when discussing in detail the metal-carbon hybridization 
levels. In brief, the behavior of the binding energies of 
the substitutional 3d transition metal comes from two 
competing effects: 

(i) as the 3d shell becomes occupied and moves to lower 
energies the hybridization with the carbon vacancy states 
near the Fermi energy (Ep) is reduced, which decreases 
the binding energy; 

(a) the transition from Mn to late transition metals is 
accompanied by a shift of the metal-carbon bond length 



of '-^0.1 A, which increases the carbon-metal interaction 
and, correspondingly, the binding energy. 

Binding energies for noble metals are considerably 
smaller than for transition metals and mirror the reverse 
behavior of the bond lengths: 3.69, 1.76 and 2.07 eV, 
respectively, for Cu, Ag and Au. The smallest binding 
energy (^1 eV) among the metals studied here is found 
for Zn, withh filled s-d electronic shells. 



C. Spin moments 

The spin moments of substitutional transition and no- 
ble metals in graphene are shown in FigJ31 Again they 
are in agreement with the results of Ref.|23- However, we 
advance a simple model to understand the observed be- 
havior which was not presented in that reference. One of 
the fundamental results of our study is a detailed model 
of the bonding and electronic structure of substitutional 
transition metals in graphene. As we will see below, the 
evolution of the spin moment can be completely under- 
stood using such model. In brief, we can distinguish dif- 
ferent regimes according to the filling of electronic lev- 
els of different (bonding, non-bonding and antibonding) 
character: 

(i) all the carbon- metal bonding levels are filled for Sc 
and Ti and, correspondingly, the spin-moment is zero; 

(ii) non-bonding 3d levels become populated for V and 
Cr giving rise to a spin moment of, respectively, 1 and 

2 /is with a strong localized d character; 

{Hi) for Mn one additional electron is added to the 
antibonding d^^ level and the spin moment increases to 

3 Mb; 

(iv) finally, for Fe and heavier atoms all the non- 
bonding 3d levels are occupied and the spin moment os- 
cillates between and 1 /.ts as the antibonding metal- 
carbon levels become occupied. 

The sudden decrease of the spin moment from 3 hb 
for Mn to /is for Fe is characterized by a transition 
from a complete spin-polarization of the non-bonding 3d 
levels to a full ocuppation of those bands. However, this 
effect depends on the ratio between the effective electron- 
electron interaction within the 3d shell and the metal- 
carbon interaction (see Sec. I VII[) . As we will see be- 
low, if the hybridization with the neighboring atoms is 
artificially reduced, for example by increasing the Fe-C 
distance, Fe impurities develop a spin moment of 2 fiB- 
Our results also show that it is also possible to switch on 
the spin moment of Fe by changing the effective electron- 
electron interaction within the 3c? shell. This can be done 
using the so-called GGA-fU method. For a large enough 
value of U (in the range 2-3 eV), Fe impurities develop a 
spin moment of 1 /is. This will be also explained in de- 
tail in Sec. lVIII For the time being we just point out that 
this behavior is unique to Fe: using similar values of U 
for other impurities does not modify their spin moments. 

At the level of the GGA calculations Fe constitutes the 
border between two different characters of the spin mo- 
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ment associated with the substitutional metal impurities 
in graphene: 3d magnetism for V-Mn and a "defective- 
carbon" -like magnetism for heavier atoms. For Co, Ni, 
the noble metals and Zn the electronic levels close to the 
Ep have a stronger contribution from the carbon nearest- 
neighbors and resemble the levels of the isolated Da/i car- 
bon vacancy. In particular, MuUiken population analysis 
show that the spin moment of the noble metals impurities 
has a dominant contribution for the three nearest carbon 
neighbors (see Table H]). For Zn two electrons occupy a 
two-fold degenerate level reminiscent of the E sp level of 
the unreconstructed carbon vacancy (see Sec. lIVp . As a 
consequence, the system suffers a Jahn- Teller distortion 
and has a zero spin moment. However, it is possible to 
stabilize a symmetric configuration (Znc3„ ) with a mo- 
ment of 2 jLts and only slightly higher in energy. 



IV. UNRECONSTRUCTED Dgfe CARBON 
VACANCY 

We have seen in the previous summary of results that, 
as substitutional impurities in graphene, most of the 
metal atoms studied here present a threefold symmetrical 
configuration. For this reason we have found particularly 
instructive to analyze their electronic structure as the re- 
sult of the hybridization between the atomic levels of the 
metal atoms with the electronic levels associated with an 
unrelaxed D^h symmetrical carbon vacancy. 

Figure. [3] (a) shows the electronic structure of such 
^3h carbon vacancy as calculated using 4x4 graphene 
supercell, while panel (b) presents a simplified scheme 
that highlights the defect levels associated with the va- 
cancy and indicates their different character and symme- 
try. The defect levels of the D3/1 vacancy can be eas- 
ily classified according to their sp or character and 
whether they transform according to A or E-type repre- 
sentations. Close to the Ej? we can find a fully symmetric 
A pz level (thus belonging to the Aj irreducible repre- 
sentation) and two degenerate (at F) defect levels with 
E symmetry and sp character (E' representation). Ap- 
proximately 4 eV below Ej? we find another defect level 
with A sp character (A'^ representation). 

It is interesting to note that it is not possible to iden- 
tify any localized defect level with E pz (E") character. 
This is due to the strong coupling with the delocalized 
states in the graphene layer and contrasts to the case of 
the A Pz level. The A pz level lies very close to E^?, where 
the density of states is low, and due to its A symmetry 
cannot appreciably couple to the delocalized Pz states of 
graphene in that energy range. On the contrary, the E 
Pz combinations present a very strong hybridization with 
the rest of the states of the graphene layer. Indeed, an in- 
spection of the Projected Density of States (PDOS) (see 
the scheme in Fig.|4](b)) reveals that the spectral weight 
associated with such E symmetry linear combinations of 
Pz orbitals of the carbon atoms surrounding the vacancy 
extend over the whole valence and conduction band of 



graphene. It is important to take the last observation 
into account when developing a model of the electronic 
structure for the metal substitutionals in graphene. Note 
that we need to have the correct number of electrons from 
carbon available for forming localized (covalent) bonds. 

The three carbon atoms around the vacancy provide 
three unpaired electrons associated with the unsaturated 
sp lobes and three electrons coming from the Pz orbitals. 
As shown in Fig. d] (b), two of these electrons stay in pz 
states delocalized over the graphene layer while the other 
four electrons fill the A sp and A pz levels localized at 
the vacancy. 



V. ANALYSIS OF THE ELECTRONIC 
STRUCTURE 

We now turn to the problem of the electronic structure 
of 3d transition and noble metal atoms as substitutional 
impurities in graphene. We first present a model of the 
hybridization between carbon and metal levels and, sub- 
sequently, we show that this model allows to understand 
in detail the band structures obtained in our calculation 
for all the metals. 



A. Sc and Ti: filling the vacancy-metal bonding 
levels 

Figure[5](a) presents a schematic representation of the 
hybridization of the 3d levels of Ti with those of an un- 
reconstructed D^h carbon vacancy. We only consider ex- 
plicitly the 3d states of the metal atom since our calcula- 
tions show that, for transition metals, the main contribu- 
tion from 4s orbitals appears well above Fip. Due to the 
symmetric position of the metal atom over the vacancy 
the system has a Cs^ symmetry and the electronic levels 
can still be classified according to the A or E irreducible 
representations of this point group. Of course, metal and 
carbon vacancy states only couple when they belong to 
the same irreducible represention. Thus, occupied A pz 
and A sp vacancy levels can only hybridize with the 3dz2 
orbitals (Ai representation), while all the other 3d metal 
orbitals can only couple to the unoccupied E sp vacancy 
levels. 

With these simple rules in mind and taking into ac- 
count the relative energy position of carbon and metal 
levels, that changes as we move along the transition metal 
series, we can propose a model of the electronic struc- 
ture of substitutional transition metals in graphene as 
represented in Fig. [5] (a) and (b) . There are three local- 
ized defect levels with Ai character and three twofold- 
degenerate levels with E character. Two of these E levels 
correspond to bonding-antibonding sp-d pairs, while the 
third one corresponds to 3d non-bonding states. For Sc- 
Mn the three Ai levels can be pictured as a low lying 
bonding level with A sp-dz^ character and a bonding- 
antibonding pair with A Pz-dz^ character. 
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FIG. 5: (Color online) (a) Scheme of the hybridization be- 
tween the 3d levels of Ti and the localized impurity levels of 
the Dsh C vacancy. Only d levels of Ti are taken into account 
since our calculations show that, at least for transition met- 
als, the main contribution from s levels appears well above 
Ef. Levels with A symmetry are represented by gray (green) 
lines, while those with E symmetry are marked with black 
lines. The region close to Ef is highlighted by a (red) square, 
(b) Schematic representation of the evolution of the electronic 
structure near Ef for several substitutional transition metals 
in graphene. The spin moment (S) is also indicated. Substi- 
tutional Sc impurities act as electron acceptors, causing the 
p-doping of the graphene layer. 



Therefore, as shown in Fig. [5] we have four metal- 
vacancy bonding levels (two A and one doubly- 
degenerate E levels) that can host up to eight electrons. 
Ti contributes with four valence electrons, and there 
are four electrons associated with the localized carbon- 
vacancy levels. Thus, for Ti the bonding states are com- 
pletely occupied. Consequently, Ti presents the highest 
binding energy among all 3d transition metals and has a 
zero spin moment. 




FIG. 6: (Color online) Calculated band structure of sub- 
stitutional Sc (a) and Ti (b) impurities in 4x4 supercell of 
graphene. Open circles indicate the contribution from 3d^2 
orbitals of the metal atom and C 2pz orbitals of the neigh- 
boring C atoms. Solid circles indicate contributions from the 



rest of the 3d orbitals and C 2px 
are referred to the Fermi energy. 



and 2py orbitals. Energies 



The situation for Sc ougth to be discussed in detail. 
As Sc has three valence electrons, in principle we could 
expect an incomplete filling of the metal-vacancy bond- 
ing levels and a spin moment of '^1 ^b- However, in our 
model the highest bonding state (with A Pz-d^'^ char- 
acter) appears below Ei? and the Sc impurity can act 
as an acceptor impurity. Our calculations show that 
this is indeed the case. The Sc-vacancy Pz-d^^ impu- 
rity level captures an electron from the extended states 
of the graphene layer which becomes p-doped. In total, 
the substituional Sc-graphcnc system docs not show any 
spin polarization. 

We can now contrast the expectations from our model 
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TABLE I: MuUiken population analysis ol the spin moment in 
the central metal impurity (Sa/) and the carbon nearest neigh- 
bors (So) for different substitutional impurities in graphene. 
Stot is the total spin moment in the supercell. 





Sm{iib) 




Stot (/x_b) 


V 


1.21 


-0.09 


1.0 


Cr 


2.53 


-0.20 


2.0 


Mn 


2.91 


-0.10 


3.0 


Co 


0.44 


0.06 


1.0 


Cu 


0.24 


-0.03, 0.31, 0.31 


1.0 


Ag 


0.06 


-0.31, 0.54, 0.54 


1.0 


Au 


0.16 


-0.28, 0.50, 0.50 


1.0 


Znc3„ 


0.23 


0.37 


2.0 



with actual calculations. Figure [6] shows the band struc- 
ture of Sc (a) and Ti (b) close to the E^. As expected, 
the main contribution from the 3d shell is found above 
Ep. Below Ep we find one defect band with p^-d^^ char- 
acter and two bands (degenerate at T) with sp-d char- 
acter. These bands are in close correspondence with the 
bonding A and E levels appearing in our model. In the 
case of Ti the E^? is located inside a gap of ^^0.5 eV 
that opens at K point in the Brillouin zone. This gap 
appears due to the relatively small 4x4 supercell used 
in these calculations and is reduced when larger super- 
cells are used. Thus, the filling of the graphene extended 
bands is not appreciably changed by substitutional dop- 
ing with Ti. For Sc the situation is different. As shown 
in Fig. [U (a) Ej? moves away from the K point, the Sc- 
vacancy complex captures one electron and the graphene 
layer becomes doped with holes. 

Regarding the unocuppied bands, the 3c? contribution 
for Sc above Ej? appears quite broadened due to the 
strong hybridization with the graphene states. Indeed, 
the defect levels are somewhat difficult to identify and to 
correlate with our model. One exception is a flat band 
with strong dz^ character appearing at ~1.5 eV that, due 
to its symmetry, does not couple so efficiently with the 
host states. The case of Ti is much easier to interpret 
in terms of the simplified model presented in Fig. [5] (a). 
In particular, we can find two bands at '^0.6 eV with 
strong dxy and d,j.2_y2 contribution that correspond with 
the non-bonding d impurity levels, and one band with 
dz'i character at ~0.8 eV corresponding with the A pz- 
dz2 antibonding level. Around 2.6 eV we can also find 
the two E sp-d antibonding defect levels, although in this 
case much more hybridized with the host. 



B. V, Cr and Mn: 3d magnetism 

As we have seen above, the metal-vacancy bonding lev- 
els are completely filled for substitutional Sc and Ti and, 
as a consequence, the spin moment associated with these 
impurities is zero. However, following our model in Fig- 




K TMK TK TMK F 

FIG. 7: (Color online) Calculated band structure of substitu- 
tional V (panels (a) and (b)), Cr ((c) and (d)) and Mn ((e) 
and (f)) impurities in a 4x4 supercell of graphene. Open and 
filled circles indicate, respectively, the contribution from 3d^2 
and the rest of the 3d orbitals of the metal and, therefore, 
also indicates levels with A and E symmetries. Energies are 
referred to the Fermi energy 
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ures[5](a) and (b), as we move along the transition metal 
series the number of valence electrons increases and the 
non-bonding 3d impurity states start to become occu- 
pied. Hence, due to the strong atomic character and lo- 
calization of these states, the system develops a non-zero 
spin moment. 

Figure [7] shows the calculated band structures for V, 
Cr and Mn impurities in a 4x4 supercell of graphene. 
The full calculations and the predictions of our simplified 
model agree remarkably, at least in the neighborhood of 
Fip. For V and Cr this correspondence is particularly 
evident: one and two electrons, respectively, occupy the 
degenerate non-bonding E d levels. These non-bonding 
levels have a dominant contribution from the Sd^y and 
3dx2-y2 orbitals of the metal atoms. As expected, the 
spin moments associated with these impurities are, re- 
spectively, 1 and 2 /is. The strongly localized character 
of these spin moments is corroborated by the MuUikcn 
population analysis shown in Table IH This analysis in- 
dicates that the spin moment is mainly localized at the 
metal impurity. The contribution from the neighboring 
carbon atoms is much smaller and has the opposite sign. 
The localized character of the moment is also consistent 
with the relatively large values of the spin splitting of the 
impurity bands. From Figure [7] we calculate a spin split- 
ting of ~0.9 eV for V and almost 2 eV for Cr for the E 
d levels at E^?. These splittings are comparable to those 
of d-electrons in magnetic bulks, in the order of 1 eV. 

Figure [7] (e) and (f ) show the majority and minority 
spin band-structures for Mn substitutionals in graphene. 
In addition to the non-bonding 3d levels, the antibond- 
ing A Pz-d^-i defect level also becomes occupied and spin- 
polarized. This level has an important contribution from 
the 3d^2 state of Mn (given by open symbols). There- 
fore, it is relatively localized within the Mn atom and 
presents a significant tendency towards spin polarization. 
Indeed, our calculations show that Mn substitutions in 
graphene give rise to a spin moment of 3 ^b- The Mul- 
liken decomposition for Mn in Table [T] confirms again the 
localized character of this spin moment. However, some 
differences respect to V and Cr are also found. In those 
cases, the moment associated with the metal atoms was 
always larger that the total moment and the only sig- 
nificant additional contributions came from the nearest 
carbon neighbors. However, for Mn the atomic moment 
is somewhat smaller (2.91) than the total moment (3.0). 
Taking into account the contribution from the nearest 
carbon neighbors (-0.3), a moment of ^0.4 /is is assigned 
to carbon atoms that are further away from the defect. 
This indicates a slightly more delocalized character of the 
spin moment of Mn, since for V and Cr the "long range" 
contribution was smaller than 0.1 hb- From the band 
structures in Fig. [7] (e) and (f) we obtain a spin split- 
ting of ~2.1 eV for the non-bonding d levels of the Mn 
impurity , similar to the case of Cr. The spin splitting 
for the antibonding A Pz-d^i state has a smaller value of 
~1.5 eV, indicative of its larger spatial extension. 

From the simple scheme presented in Fig. [5] (a) we can- 



not completely determine the value of the spin moment 
of the Mn impurity. It can be 3 iib as found in our 
first principles calculations and schematically depicted 
in Fig. [5] (b). However, a magnetic moment of 1 /is is 
also a possible answer. In the latter case, the additional 
electron in Mn with respect to Cr could populate one of 
the minority-spin non-bonding d impurity levels instead 
of the antibonding A Pz-dz2 level. In such situation, the 
spin moment is determined by a delicate balance between 
the onsite exchange energy within the 3d shell, and the 
energy cost (Ae^^) to promote one electron from the non- 
bonding d levels to the higher energy antibonding A p^- 
d^2 state. Note that the electron-electron repulsion is 
also reduced when the electron moves into the less local- 
ized A Pz-dz2 level. An estimate of the exchange energy 
can be obtained from the spin splitting (A5) of the de- 
fect levels nearby E^. The relative position of the 3d 
states respect to the A pz level of the carbon vacancy 
and the interaction matrix element between these levels 
determine Ae^d- Thus, within our GGA-DFT calcula- 
tions we can expect the high spin solution to be favored 
approximately when A5 > AeAd- From the band struc- 
tures in Fig. [7] (e) and (f) we obtain AcAd 1-0 eV, 
which is smaller than the values of Ag discussed pre- 
viously and, therefore, is consistent with the calculated 
moment of 3 I-Ib for Mn. Fe impurities considered below 
in detail present a similar situation where two spin con- 
figurations are possible. However, in the case of Fc the 
low spin (spin compensated) solution is preferred at the 
level of DFT-GGA calculations as a result of the stronger 
metal-carbon hybridization. 

In short, both the results of the calculations and the 
expectations based on our model of the metal-vacancy 
bonding point towards a very strong 3d character of the 
defect levels appearing nearby of Ei;- for V, Cr and Mn 
substitutional impurities in graphene. The filling of these 
localized levels favors high spin solutions in accordance 
with the first Hund's rule of atomic physics. Thus, we 
can picture the appeareance of spin polarization for V, Cr 
and Mn substitutionals in graphene as "standard" d-shell 
magnetism. 



C. Fe, Co, Ni: strong contribution from the carbon 
vacancy levels 

As we have seen in the previous section the defect lev- 
els appearing in the neighborhood of E^ associated with 
the presence of V, Cr and Mn substitutional impurities 
in graphene have a strong 3d character. Consequently, 
these impurities exhibit large spin moments. However, 
when increasing the atomic number along the transition 
series, the atomic 3d levels move to lower energies and 
we enter a different regime: the defect states nearby 
have a predominant contribution from the carbon atoms 
neighboring the metal impurity. For 3d late transition 
metals heavier than Fe, the noble metals and Zn we can 
stablish a strong link between the electronic structure 
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(a) 




FIG. 8: (Color online) (a) Similar to Fig. [5] (a) for the case of 
Ni. The region close to the Fermi level is indicated by a (red) 
circle, (b) Scheme of the levels close the Ef for Fe, Co, Ni, 
noble metals and Zn. For Fe, in addition to the antibonding 
metal-vacancy levels, we have also included the non-bonding 
d levels that also appear quite close to Ef. The noble metals 
slightly break the Cs^ symmetry. For Zn there are two so- 
lutions: a high spin solution that preserves symmetry and a 
distorted one with zero spin moment. 



of these impurities around E^? and that of the unrecon- 
structed D^h carbon vacancy. 

A detailed scheme of the hybridization of the Ni 3d 
states with those of the carbon vacancy, and the result- 
ing electronic structure, is presented in Fig. [5] (a). All the 
bonding and non-bonding metal-vacancy states are filled 



in this case and the levels appearing closer to E^? have an- 
tibonding character and a strong contribution from the 
three carbon nearest-neightbors. Closely below E^? we 
find a level with A character and a strong C 2pz contribu- 
tion and a smaller weight in the Ni 3dz2 state. Above Ep 
there are two degenerate levels with E character mainly 
coming from the C 2sp lobes and hybridized with the Ni 
3dxz and 3dyz orbitals. The resulting electronic struc- 
ture strongly resembles that of the isolated T>3h (unre- 
constructed) carbon vacancy as can be checked by com- 
paring the band structures in Fig. [4] (a) and Fig. [9] (d). 
The main difference stems from the slightly higher po- 
sition of the unoccupied E levels in the case of the Ni 
impurity. This upward shift is due to the antibonding 
interaction with the d states of Ni and contributes to the 
stability of the C^y spin-compensated solution. In the 
case of the D^h carbon vacancy, the E C 2sp levels lie 
closer to Ep and make the system unstable against spin 
and structural distortions. This clear connection between 
the electronic structure of Ni and that of the Da/; van- 
cancy was already emphasized in Ref. [22I . 

Thus, our GGA calculations predict Ni substitutional 
to be a closed shell system with zero spin moment. One 
could expect that a better description of the electron- 
electron interaction within the 3d shell would enhance 
any tendency of the system towards a magnetic instabil- 
ity. For this reason we have performed GGA-I-U calcula- 
tions (using the Vasp code) with values of U up to 4.5 eV. 
However, Ni substitutional impurities in flat graphene al- 
ways remain non- magnetic in the calculations. 

In Ref. m we proposed a different way to switch on 
the magnetism of Ni substitutionals: the possibility of 
spin-moments induced by curvature. The idea consists in 
lifting the degeneracy of the two unoccupied E sp-d levels 
by applying a structural distortion. If the distortion is 
large enough one of these levels close toEp becomes par- 
tially populated and, due to its small band width, spin 
polarized. We have checked that the curvature of the 
carbon layer in (n,n) nanotubes with n ranging between 
4 and 8 induces spin moments as large as 0.8 /is per Ni 
substitutional impurity. The spin moment in these sub- 
stitutionally Ni-doped nanotubes is strongly dependent, 
not only on the layer curvature, but on the density and 
arrangement of defects in the tube. We have recently 
demonstrated that a similar switching of the magnetic 
moment can be obtained in flat Ni-doped graphene by 
applying adequate structural distortions, thus providing 
a simple way to control the spin of this system4i 

As suggested by our scheme in Fig. |S] (b), we can now 
try to understand the electronic structure of Co and Fe 
impurities from that of the Ni substitutional but remov- 
ing, repectively, one and two electrons. According to this 
image Co substitutionals should present a spin moment 
of 1 fiB^ while Fe substitutionals should be non- magnetic. 
This is indeed confirmed by our GGA calculations. 

Figures [5] (b) and (c) show, respectively, the majority 
and minority spin band structures of the Co impurity. 
In the neighborhood of E^? we find a spin-polarized band 
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FIG. 9: (Color online) Like Fig.[Z]for Fe (a), Co (b) and (c), Ni (d), Cu (e) and (f), Zn (g) and symmetric Zncg^ (h) and (i) 
substitutional impurities in graphene. 



associated with the antibonding A Pz-dz^ impurity level. 
The spin splitting of this band is ^0.5 eV. The hybridiza- 
tion character of this level is confirmed by the Mulliken 
analysis in Table. [H Only a contribution of 0.44 /is to 
the total spin moment comes from the Co atom. The 



relatively delocalized character of the A Pz-dz^ level also 
becomes evident. Only a moment of 0.18 fiB comes from 
the three carbon nearest neighbors, while 0.38 fiB comes 
from carbon atoms at larger distances. The slow distance 
decay of the A Pz-dz^ defect level translates into quite 
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TABLE II: Structural parameters for substitutional noble 
metals in graphene. 





dc-M(A) 


h.(A) 


en 


Cu 


1.93, 1.90, 1.90 


1.40 


88.9, 88.9, 95.2 


Ag 


2.23, 2.19, 2.19 


1.84 


71.7, 71.7, 76.7 


Au 


2.09, 2.12, 2.12 


1.71 


78.0, 78.0, 81.6 



strong and long-range magnetic interactions between mo- 
ments associated with neighboring Co defects. Indeed, 
the peculiar electronic structure of the Co impurities has 
important consequences for the magnetism of this sys- 
tem: couplings show a complex dependence with distance 
and direction, while the total spin-moment is determined 
by the number of Co substitutions in each sublattice of 
the graphene layer. We refer the interested reader to 
Ref.ii. 

Figure [9] (a) presents the GGA band structure for a Fe 
substitutional defect in a 4x4 supercell of graphene. Sim- 
ilar results are found using a larger 8x8 supercell. This 
band structure is again in reasonable agreement with the 
simple model presented in Fig.[8l The non-bonding d lev- 
els are completely filled and appear ^0.4 eV below Ep in 
the vicinity of F. The A Pz-d^^ level is mostly unoccupied 
and close to 0.2 eV above Ep near F. However, we can 
see that nearby the K point the Pz-d^^ band becomes par- 
tially occupied, indicating a small charge transfer from 
the dispersive tt bands of graphene to the defect. Mul- 
liken analysis also reflects a small charge accumulation of 
^^0.16 electrons in Fe. In spite of this small partial popu- 
lation, the spin compensated solution is the most stable 
for Fe substitutionals at the GGA level. 

As already pointed out in Sec. IIII CI the magnetic be- 
havior of the Fe imputiy is a consequence of a delicate 
balance between the onsite electron-electron interaction 
and the metal-carbon hybridization. For this reason, we 
devote a whole Section below (Sec. IVIip to explore how 
the band structure and spin moment of Fe substitution- 
als are modified when these factors are independently 
controlled by changing the Fe-C bond length and using 
the GGA-f U approximation to describe the effects of the 
electron-electron repulsion within the d-shell. 



D. Noble metals 

In the previous sections we have seen that some of the 
traditional ferromagnets, like Fe and Ni, become non- 
magnetic as substitutional impurities in graphene. Also 
quite surprisingly, we find here that substitutional impu- 
rities of noble metals are magnetic with a spin moment 
of 1 /is . The reason for this behavior becomes clear once 
the electronic structure of these defects is understood. 

The band structures can be found in Fig. [9] (e) and 
(f) for Cu and in Fig. [10] for Ag and Au. In agreement 
with the predictions of our simplified model in Fig [8l we 
can see that in the case of the noble-metal impurities 




FIG. 10: (Color online) Like Fig. [7] for Ag and Au. The 
similarities with the band structure of Cu in Fig.|9]are evident. 



the two-fold degenerate E sp-d antibonding levels^i are 
now occupied with one electron. The system undergoes a 
small structural distortion that removes the degeneracy 
of these levels and the unpaired electron becomes spin po- 
larized. Therefore, substitutional impurities of the noble 
metals in graphene exhibit a spin moment of 1 ■ Rela- 
tivistic effects are known to be much more important for 
Au than for Ag and Cu. Although we have not included 
spin-orbit coupling in our calculations, scalar relativistic 
effects are taken into account in the construction of the 
pseudopotentials. However, the similarities between the 
electronic structure of all three noble metals are evident, 
which indicates that bonding and magnetic behavior is 
mainly dictated by the number of valence electrons. 

Structural parameters for all noble metals impurities 
can be found in Table [TTl For Cu and Ag one of the metal- 
carbon bond lengths is slightly larger than the other two. 
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TABLE III: Structural parameters and binding energies for 
substitutional Zn impurities in graphene for the symmetric 
Csi, and most stable distorted configurations. 





Ei3 (eV) 


dc-M (A) 


(A) 


en 


Zn 


0.91 
1.07 


1.99 
2.06, 1.89, 1,89 


1.67 
1.54 


87.9 
88.3 88.3 103.9 



whereas for Au one is shorter than the others. However, 
the distortions are rather small with variations of the 
bond lengths smaller than 2 %. The differences intro- 
duced by the larger scalar relativistic effects of Au mainly 
reflect in the slightly smaller metal-carbon bond length 
for this metal as compared with Ag. 

Table U shows the distribution of the spin moment 
among the metal atom and the nearest carbon neighbors. 
We can see that the contribution from the metal atom is 
almost negligible, particularly in the case of Ag. This 
can be expected from the lower energy position of the d 
shell in the case of the noble metals as compared with 
transition metals. Although slightly hybridized with the 
p shell of the metal impurity, the defect states nearby E p 
in this case are mainly coming from the carbon neighbors. 
Still the spin moment is rather localized in the complex 
formed by the metal atom and its three nearest neigh- 
bors, which for Cu and Ag contributes with a moment 
of 0.83 /is, and up to 0.88 hb for Au. The contribu- 
tion from the rest of the graphene layer is much smaller 
than, for example, in the case of Co. This reflects the 
dominant contribution from the relatively localized car- 
bon sp lobes in the description of the defect states nearby 
Ei? for these impurities. This analysis reinforces the link 
with the electronic structure of the unreconstracted D^h 
carbon vacancy in graphene as presented in Sec. IV CI In 
fact, we can picture the main role of late transition- and 
noble-metals substitutionals in graphene as to stabilize 
the structure of the carbon monovacancy, that otherwise 
will severely reconstruct, and to change its charge state. 



VI. JAHN-TELLER DISTORTION OF 
SUBSTITUTIONAL ZN 

For Zn impurities a second electron is added to the two- 
fold degenerate E sp-d shell. Under these circumstances 
two scenarios are possible: (i) a non-magnetic solution in 
which the system has undergone a Janh- Teller-like dis- 
tortion, or [a) a high-spin solution that maintains the 
symmetric C3.U geometry of the defect. The relative en- 
ergy of both solutions depends on the balance between 
the energy gain associated with the distortion and the ex- 
change energy of the electrons. Both types of solutions 
are obtained in our density functional calculations for Zn 
substitutional impurities in graphene. 

The details of the structure and the binding energies 
of the Zn impurity are presented in Table IIIII as calcu- 
lated with Siesta. Very similar results are obtained for 



both configurations using Vasp. The distorted config- 
uration presents one larger Zn-C bond (by ~3.5%) and 
two shorter bonds (^5%) compared with the bond length 
(1.99 A) of the undistorted geometry. The distorted con- 
figuration is more stable by 160 meV (120 meV using 
VASP). This rather small energy difference between the 
two configurations might point to the appearance of non- 
adiabatic electronic effects at room temperature. 

The band structure for both configurations of Zn sub- 
stitutionals can be found in Fig. [51 Again they con- 
firm the model presented in Fig. [51 The distorted Zn 
[Fig. [51 (g)] breaks the degeneracy of the E sp-d levels: 
one of them appears fully occupied '^0.8 eV below Ep, 
while the other appears a fews tenths of eV above Ep. 
For the C^y Zn impurity both E sp-d bands are degen- 
erate and the splitting between majority and minority 
levels is ^0.71 eV. Table U shows the MuUiken popula- 
tion analysis of the spin moment for the C^y Zn systems. 
As in the case of the noble metals the contribution from 
the three nearest carbon neighbors is the most important. 
However, the contribution of Zn is somewhat larger and 
the total spin moment is more delocalized with a contri- 
bution form the rest of the graphene layer of ~0.66 fiB- 
We should note that here we have one additional electron 
as compared to the noble metal systems. 



VII. FE SUBSTITUTIONALS: COMPETITION 
BETWEEN INTRA-ATOMIC INTERACCIONS 
AND METAL-CARBON HYBRIDIZATION 

We have already pointed out that Fe substitutionals in 
graphene occupy a rather special place at the border be- 
tween two well defined regimes (see Fig.[3|): (i) the strong 
3d character of the defect levels nearby Ep and large spin 
moments found for V, Cr and Mn impurities, and (ii) the 
larger carbon character of those electronic levels and the 
small oscillatory spin-moments of Co, Ni and the noble 
metals. GGA calculations [Fig. [51 (a)] locate Fe impuri- 
ties within the second group, with all the 3d non-bonding 
levels fully occupied and non-magnetic. Thus, the spin 
moment drops from 3 fiB for Mn impurities to zero for 
Fe, showing a quite discontinuous behavior as a function 
of the number of valence electrons. In analogy with the 
standard Slater-Pauling rule for transition metals, one 
could expect to find a more gradual decrease of the spin 
moment as the number of valence electrons is increased, 
i.e. Fe would have a moment of 2 fiB- In the present 
section we study in detail why the non-magnetic solution 
is more stable for Fe. 



A. Key parameters: metal-carbon hopping and 
intra-atomic Coulomb interactions 

Figure [TT] shows a scheme of the electronic structure of 
both Mn and Fe defects. Depending on how electrons are 
arranged among the A Pz-d^^ antibonding level and the 
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FIG. 11: (Color online) Scheme of the electronic levels near 
Ef for Mn and Fe substitutionals in graphene as deduced 
from our model of bonding and GGA calculations. Aexd is 
the energy cost to promote an electron to the antibonding 
A p2-d^2 hybridization level from the non-bonding 3d states. 
The magnitude of AeAd, relative to that of the spin-splitting 
As of these defect levels, is crucial to determine the spin state 
of these impurities. 



non-bonding 3d states, Mn can exhibit spin moments of 
3 hb OT 1 fj,B, while Fe can have a moment of 2 /^s or can 
be non-magnetic. At the GGA level Mn prefers the high- 
spin configuration, while the low-spin one is more stable 
for Fe. As commented in Section IV Bl the relative sta- 
bility of the different spin states depends on the balance 
between the effects of Coulomb repulsion and exchange 
within the 3d levels and the relative energy position of 
the impurity levels given by AeAd (see Fig. [TT|). The 
hybridization with the neighboring C atoms is crucial in 
this interplay since it influences (z) the degree of local- 
ization of the defect levels and the screening of Coulomb 
interactions, which modify the spin splitting of the elec- 
tronic levels A5, and (ii) the value of AcAd through the 
effective hopping parameter between the 3c? states and 
the A pz carbon vacancy level. Next, we shall deal with 
both aspects for subtitutional Fc in graphene. 




dc-Fe ( A) 



^b^c-Fe=2-08 A up 





rK r M K 



FIG. 12: (Color online) Panel (a) shows the spin moment 
of a substitutional Fe impurity in graphene as a function of 
the C-Fe bond length. Equilibrium position corresponds to 
dc-Fo=1.78 A. Panels (b) and (c) show the band structure 
for dc-Fo=2.08 A. Symbol code similar to Fig. [T] Energies 
referred to the Fermi level (indicated by a dashed line). 



1. Changing Fe-graphene hopping with distance 

By tuning the interaction with the host structure it 
should be possible to change the spin moment of these 
impurities. We can modify the hopping by artificially 
changing Fe-graphene distance. The results of such cal- 
culation are shown in Fig.[T3(a). While maintaining the 
three-fold symmetry of the system, we have performed 
a series of calculations by progressively increasing the 
height of the Fe atom over the graphene layer. Increas- 
ing the C-Fe bond length by ~9% we observe an abrupt 



jump of the spin moment from zero to 1.5 fiB- The spin 
moment continues to rise and saturates at a value of 
2.0 fiB for dc-Fe'^2.07 A. This convincingly shows that 
the metal-carbon hybridization is the key parameter to 
explain the non-magnetic state of Fe substitutional im- 
purity in graphene. When increasing the C-Fe distance 
we mainly decrease AeAd- Thus, we reduce the energy 
penalty for promoting electrons from the non-bonding d 
to the A Pz-dz2 defect levels. At the same time we also 
increase the atomic character of the non-bonding d states 
and reduce the effect of the screening due to the electrons 
in the graphene layer. This AeAd reduction promotes the 
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FIG. 13: (Color online) Projected density of states (PDOS) 
and total density of states (DOS) for a Fe substitutional im- 
purity in graphene calculated with VASP using GGA and 
GGA+U with \J—2.4 eV. Positive (negative) values for ma- 
jority (minority) spin. When U=0 the system shows zero 
spin-polarization and the Fe 3d^2 state is partially occupied. 
A large enough Hubbard U forces an integer population of the 
Fe 3d^2 state and produces a spin moment of 1 /is. 



electron-electron repulsion within the 3d states of Fe and, 
eventually, stabilizes the solution with 2.0 ^b- 

Figures [T^] (b) and (c) present the band structure for 
the high-spin state of Fe, where the C-Fe bond length has 
been elongated up to 2.08 A. We notice the differences 
with the spin-compensated ground state in Fig. [3] (a). 
On the one hand, one electron is promoted from the 
non-bonding d states to the more delocalized antibond- 
ing A Pz-dz2 impurity level in order to reduce the effect 
of electron-electron repulsion. On the other hand, the 
spin moment is maximized in accordance with Hund's 
first rule. 

Therefore, we have seen that the non-magnetic char- 
acter of the Fe substitutionals in the GGA calculations 



is due to the larger interaction with the graphene layer, 
as compared for example with the Mn impurity. This is 
consistent with the fact that Fe, together with Co and Ni, 
presents the smaller carbon metal bond length among the 
whole series of 3d transition metals. Fe impurities also 
have one of the largest binding energies. By artificially 
reducing this interaction it could be possible to obtain 
magnetic Fe substitutional impurities in graphene. 



2. Changing intra-atomic Coulomb interaction U 

Another route to explore would be to increase the size 
of intra-atomic electron-electron interactions. We have 
done that by using the so-called GGA-I-U methodology 
in which a Hubbard term is explicitly added to the DFT 
Hamiltonian and solved within the mean-field approxi- 
mation. Our results indicate that with a reasonable value 
of U (~2 eV or larger) we obtain a magnetic solution for 
Fe substitutionals. However, contrary to our initial ex- 
pectation this solution does not correspond to the 2 /zb 
high-spin solution discussed above, but to a new solution 
with 1 fiB- The key reason to understand this behavior 
is the partial occupation of the Fe 3^22 state at the level 
of GGA calculations. The state is strongly coupled 
to the delocalized A pz vacancy level appearing nearby 
and both, this hybridization and the population of 
the atomic orbital, arc strongly modified when the value 
of U is increased. 

Figure [T3l shows the results of the electronic structure 
of Fe impurities calculated using GGA+U with U=2.4 eV 
and compared with those with U=0. The two upper pan- 
els show the projected density of states (PDOS) onto 
the dxy orbital (the projection onto the d^^^y^ orbital 
is identical by symmetry), the two middle panels show 
the PDOS onto the d^^ orbital of Fe and the two lower 
panels the total density of states. At the level of GGA, 
with U=0, there is a very well defined peak closely be- 
low Ei? in the d^y PDOS corresponding to the position 
of the non-bonding d states. The spectral weight coming 
from the dy2 orbital is more spread: presents some broad 
structure around —1.5 eV corresponding to the bonding 
A Pz-dz2 defect level interacting with the valence band 
of graphene, and a narrowers peak near E^^^ with origin 
in the slightly occupied antibonding A Pz-d^-^ impurity 
state. Therefore, from the Fe 3c? states appearing nearby 
Ei? we can picture the d^y and dx^_y2 states as fully oc- 
cupied and the dz2 as partially occupied due to the larger 
interaction with the carbon neighbors. 

For a sufficiently large value of U the C 2pz-'idz2 
hybridization is overcome by the tendency of the elec- 
trons to have an integer population within the local- 
ized Zd shell. Thus, the 3^22 localizes one electron and, 
as a consequence, the system develops a 1 /is spin- 
polarization. This changes can be appreciated in Fig. 1131 
The majority-spin ?idz2 PDOS shifts downwards, while a 
strong unoccupied peak appeas around 1.5 eV in the Zdz'^ 
minority spin PDOS. Simultaneously, a more delocalized 
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level mostly coining from the 2pz orbitals of the neighbor- 
ing C atoms appears half- filled at Ej?. Note that other 
levels, such as the dxy, do not suffer such strong modifi- 
cations since they are already almost fully occupied with 
U=0. It is also interesting to note that, due to symmetry 
considerations, they couple with the graphcne layer very 
differently from the d^^ level. 



B. Revelance for recent experiments of Fe 
implantation in graphite 

There are recent reports on the paramagnetism of 
iron-implanted graphite that indicate the existence of 
local magnetic moments associated with the implanted 
Fe atoms.— This would be in contradiction with the 
present GGA results if we assume that the Fe atoms 
are incorporated to the graphene layer as substitution- 
als. However, the final geometry of the implanted atoms 
in these experiments is not known. Furthermore, a con- 
siderable amount of defects is created during the implan- 
tation process. Although it has been argued that a large 
part of the damage is healed by vacancy-interstitial re- 
combination^ the influence of these defects, specially 
interstitials and big voids, on the observed magnetic re- 
sponse can be determinant. 

Therefore, it is not fully clear if we can compare our 
calculations for Fe substitutional impurities in an other- 
wise perfect graphene layer with these experimental data. 
However, as we have shown in detail in this section, Fe 
substitutionals are very close to a transition and, depend- 
ing on the details of the calculations, it is possible to 
obtain a magnetic ground state. In particular, Fe sub- 
stitutionals develop a spin moment of 1 /is at the level 
of GGA-I-U calculations for reasonable values of the U 
parameter. This might be an indication that the non- 
magnetic ground state found in GGA calculations is a 
consequence of the limitations of the used functionals. 



VIII. CONCLUSIONS 

We have presented a DFT study of the structure, ener- 
getics, and electronic and magnetic properties of several 
metal atoms as substitutional impurities in graphene, i.e., 
bound to a carbon monovacancy in the layer. We have 
considered the cases of all 3d transition metals, noble 
metals and Zn. We have paid special attention to their 
electronic and magnetic properties and develop a simple 
model to understand the observed trends. Our model 
is based on the hybridization of the states of the metal 
atoms with those of an unreconstructed carbon vacancy. 
The main ingredients of the model are the assumption, 
after our calculations, of a three-fold symmetric bond- 
ing configuration and the approximate knowledge of the 
relative energy positions of the levels of the carbon mono- 
vacancy and the d shell of the metal impurity as we move 
along the transition series. With this model we can un- 



derstand the observed the variations of the electronic 
structure of the defect, the size and localization of the 
spin moment, and the binding energy. We have iden- 
tified three different regimes corresponding to filling of 
carbon-metal hybridization shells with different charac- 
ter: bonding, non-bonding and antibonding. 
In more detail: 

{i ) Most substitutional metal impurities present an 
almost perfecly symmetric three-fold configuration with 
Cs^ symmetry. Noble metals slightly depart from this 
perfect configuration. Only Zn presents a considerable 
structural distortion. 

[ii) For Sc and Ti the metal-carbon bonding shell is 
completely filled. Therefore, these impurities present the 
highest binding energies and zero spin moments. Sc sub- 
stitutionals act as p-dopants for graphene: each Sc im- 
purity localizes one extra electron from the carbon layer. 

[in) The non-bonding d shell becomes partially popu- 
lated for V, Cr and Mn, which develop a very localized 
spin-moment of 1, 2 and 3 /is, respectively. The binding 
energy decreases slightly as the d shell moves to lower 
energies, thus reducing its hybridization with the higher 
carbon vacancy levels. 

[iv) For Co, Ni, the noble metals and Zn, the metal- 
carbon levels are progressively populated. This gives rise 
to a oscillatory behavior of the spin moment between 
and 1 /is. The spin moments are more delocalized than 
those found for V, Cr and Mn and present a considerable 
contribution from the carbon atoms around the impurity. 
The binding energy presents a local maximum for Co, but 
suddenly drops for the noble metals and has its minimum 
for Zn. 

[v) The electronic structure, nearby E^^-, of substitu- 
tional impurities of Co, Ni, the noble metals and Zn has 
a strong resemblance to that of the unreconstructed T)zh 
carbon monovacancy We can draw an analogy between 
the electronic structure of these impurties and that of the 
unreconstructed D^h carbon monovacancy with different 
number of electrons (charge states). The spin moment of 
these impurities can be fully understand exploiting this 
equivalence. In particular, the result that noble metals 
develop a spin moment of 1 hb emerges naturally within 
this picture. 

[vi) For the Co impurity, the equivalence can be pushed 
a step further and we can draw an analogy with the 
electronic structure of a 7r-vacancy in a simple 7r-tight- 
binding description of graphene. This can be used to 
explain the peculiar behavior found for the magnetic cou- 
plings between Co substitutionals in graphene.^- 

{vii) Fe impurities occupy a distinct position at the 
border between two different regimes. Their magnetic be- 
havior stems from the competition between the carbon- 
metal hybridization and the electron-electron interaction 
within the "id shell. As a result, although Fe impurities 
are non-magnetic at the GGA level, GGA-t-U calcula- 
tions with moderate values of U (above ~2 eV) produce 
a spin- moment of 1 /is. 

[viii) We have found that the unexpected result that 
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Au substitutionals^^''^'' present a spin moment oil fin also 
holds for all noble metals, i.e., also Ag and Cu present a 
1 /is spin moment. 

{ix) We have found that the ground state of the Zn 
substitutional is non-magnetic due to a Jahn- Teller dis- 
tortion. Yet, it is possible to stabilize a symmetric con- 
figuration with a spin moment of 2fiB with a very small 
energy penalty of ~150 meV. 

In summary, substitutional impurities of metals in 
graphene present some interesting magnetic and elec- 
tronic properties and, therefore, can provide an interest- 
ing route to add functionalities or to tune the response 
of devices based on grapheme materials. Furthermore, 
recent experiments by Rodriguez-Manzo and Banharli^^ 
have demonstrated the possibility to create individual 
vacancies at desired locations in carbon nanotubes us- 
ing electron beams. This ability, in combination with 
the high stability of substitutional impurities, can open 
a route to fabricate ordered arrays of these impurities at 
predefined locations. Such devices would allow, among 
other applications, the experimental verification of the 
theoretical predictions of unsual magnetic interactions 
mediated by graphenc j^^'^^i^^ 

APPENDIX A: PSEUDOPENTENTIALS AND 
BASIS ORBITALS RADII 



the Spanish Ministerio de Educacion y Ciencia (Grant 
No. FIS2007-66711-C02-02) and the ETORTEK pro- 

TABLE IV: Cutoff radii (in Bohrs) used for the generation 
of the Troullier-Martina^ norm-conserving pseudopotentials 
used in the Siesta calculations, ri stands for the cut-off radius 
used for the / channel, while Vcore is the matching radius used 
for the construction of pseudocores in order to include non- 
linear core corrections for exchange and correlation^. 

Valence r,s (ap) (gp) rg (gp) rcore 



Sc 




3.30 


3.30 


1.44 


0.89 


Ti 


[Ar]4si3d^ 


2.96 


2.96 


1.45 


0.72 


V 




2.40 


2.79 


1.46 


0.69 


Cr 




2.51 


2.80 


1.46 


0.65 


Mn 


[Ai]-is^3d'^ 


2.51 


2.77 


1.45 


0.60 


Fe 


[Ar]4si3d'^ 


2.10 


2.10 


1.68 


0.67 


Co 


[Ar]4s2 3d^ 


2.37 


2.48 


1.68 


0.67 


Ni 


[Ar]4s''3d* 


1.85 


1.95 


1.45 


0.53 


Cu 


[Ar]4s^3d'-° 


2.33 


2.30 


1.79 


0.53 


Ag 


[Kr]5s^4d^° 


2.45 


2.58 


2.00 


0.83 


Au 


[Xe,4f*]6s^5d^° 


2.55 


2.68 


2.20 


0.93 


Zn 


[Ar]4s23dio 


2.04 


2.21 


1.66 


0.49 


C 


[He] 25^2^2 


1.25 


1.25 


1.25 





For the Siesta calculations we have used TrouUier- 
Martins norm-conserving pseudopotentials^ generated 
using the pseudization radii shown in Table IIVI of this 
Appendix. The pseudopotentials for the metal atoms in- 
clude nonlinear core corrections^^ for exchange and corre- 
lation. The pseudocore radii {vcore) have been optimized 
for each element. These pseudopententials have been 
tested to reproduce the correct results in bulk phases. 

For most elements, the shape and radii of the basis sets 
of numerical atomic orbitals (NAOs) used with Siesta 
were determined automatically by the code using an en- 
ergy shift parameter of 50 meV."^^'"*'' However, for some 
metal atoms the binding energies where slightly overes- 
timated using these basis sets, mainly due to the con- 
finement of the free atom. For those atoms we enlarged 
the radii of the basis orbitals (using smaller values of the 
energy shift parameter) until binding energies were con- 
verged within a few tens of meV. The largest radius for 
each element and I channel is shown in Table |Vl Notice 
that double-C polarized (DZP) basis of these metals in- 
clude a p shell. In this case the radii of p orbitals is taken 
equal to that of the s shell. 



TABLE V: Cutoff radii (in Bohrs) of the numerical atomic 
orbitals (NAOs) used in Siesta as a basis set. Radii used for 
the s and polarization p shell are equal. 



„NAO 



(a°) 



Co 
Ni 
Cu 
Ag 
Au 
Zn 



8.00 
10.94 
8.87 
10.48 
8.63 
9.24 



NAO 
^ d 



4.73 
6.81 
5.52 
6.52 
6.08 
5.33 
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